Modulated Scale-free Network in the Euclidean Space 
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A random network is grown by introducing at unit rate randomly selected nodes on the Euclidean 
space. A node is randomly connected to its i-th predecessor of degree ki with a directed link of 
length I using a probability proportional to kil a . Our numerical study indicates that the network is 
Scale-free for all values of a > a c and the degree distribution decays stretched exponentially for the 



other values of a. The link length distribution follows a power law: 
exactly for the whole range of values of a. 

PACS numbers: 05.10.-a, 05.40.-a, 05.50.+q, 87.18. Sn 
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Statistical properties of many different networks are 
being studied recently with much interests. Examples 
include the World-wide web (WWW) 0, the internet 
structure 0, neural networks ||, collaboration network 
etc. Broadly these networks are classified into four dif- 
ferent groups, namely, (i) Networks on regular lattices (ii) 
Random networks (^] (iii) Small-world networks (SWN) 
§ and (iv) Scale-free networks (SFN) §, §. 

It has been observed that the degree distributions of 
nodes for two very important networks e.g., World-wide 
web |l| which is a network of webpages (nodes) and the 
hyperlinks (links) among various pages and the inter- 
net network Q of routers or autonomous systems follow 
power law as: 



P(k) ~ k~ 



(1) 



These networks are called Scale-free networks and the 
exponent 7 varies between 2 and 3 for these networks. 
Barabasi and Albert (BA) proposed a simple model for 
an evolving SFN that has the following two essential in- 
gredients, namely: (i) A network grows from an initial 
set of m nodes with m < m links among them. Fur- 
ther, at every time step a new node is introduced and 
is randomly connected to m previous nodes, (ii) Any of 
these m links of the new node introduced at time t con- 
nects a previous node i with an attachment probability 
TTi (t) which is linearly proportional to the degree ki (t) of 
the i-th node at time t 



BA 



(t) ~ ki(t). 



(2) 



For BA model 7 = 3 g. 

The physical distance, or the Euclidean distance be- 
tween the nodes plays an important role in cases like 
electrical networks, internet or even in postal and trans- 
ports networks etc. In these networks one tries to mini- 
mize the lengths of the connections e.g., electrical wires, 
ethernet cables or say travel distances of postal carri- 
ers. Static networks, in which connection probabilities 
depend on the Euclidean distance have already been con- 
sidered in the context of small- world properties |J. Re- 
cently a number of interesting works in the context of 
SFN have been published 0, pj. 



Study of internet's topological structure is important 
for designing efficient routing protocols and modelling 
internet traffic. Waxman model describes internet with 
exponentially decaying link length distribution: D(£) ~ 
exp(—£/£ ) Q. Faloutsos et. al. observed the scale-free 
degree distribution of the internet S. Yook et. al. ob- 
served nodes of the router level network maps of North 
America are distributed on a fractal set and the link 
length distribution is inversely proportional to the link 
lengths . They also argued that a competition exists 
beteen the preferential attachment of the nodes and the 
weightage of the link lengths Jll|] . 

In this paper our aim is to study how a scale-free net- 
work defined on the Euclidean space behaves when the 
usual BA attachment probability as in eqn. (2) is mod- 
ulated by a link length £ dependent factor l a . Our im- 
portant observation is even for the uniform random dis- 
tribution of nodes we obtain the power law variation of 
the link length distribution for all values of a including 
the empirically observed inverse variation where as the 
Waxman's exponential behaviour is never observed. We 
argue that for a country with homogeneously distributed 
router density, our results seem to be important. 

Specifically in two-dimensions, we consider an unit 
square area on the x — y plane. Randomly selected points 
within this area are the nodes of the network. The net- 
work grows by systematically introducing one node at a 
time with randomly chosen coordinates (a;, y); < x, y < 
1 with uniform probabilities. The attachment probabil- 
ity that the new node introduced at time t would be 
connected to its i-th predecessor (0 < i < t — 1) is: 



71* (f) ~h{t)t 



(3) 



where, I is the Euclidean distance between the i-th and 
the i-th node and a is a continuously varying parameter. 

The case with a—0 is the usual BA model. For the 
negative values of a, the largest value of the modulation 
factor l a corresponds to the smallest value of I. There- 
fore, in the limit of a — > —00, only the smallest value of 
I corresponding to the nearest node will contribute with 
probability one. Similarly, for a > large I values will 
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a = a = -1 




a = -2 a = -oo 

FIG. 1: Modulated scale free networks within an unit square 
for different values of the modulation parameter a for the 
same distribution of 512 nodes. For a = +oo and — oo a newly 
introduced node is linked only to its farthest and nearest pre- 
decessors respectively, where as for a = it is connected to 
one of the previous node according to the BA rule. 

be more probable and the limit of a — > +oo corresponds 
to only non-zero contribution from the furthest node. 

We start with only one initial node corresponding to 
m = m — 1 and connect a new node to only one of its 
previous nodes. Therefore, every node has only one out- 
going link (k OU i = 1) but can have a number of incoming 
links (ki n ). The network thus formed has a tree struc- 
ture, without any loops. Similar to BA model we expect 
that the main results of our model should be robust with 
respect to the value of m used [|| . 

The link lengths of this network vary over a wide range. 
We define D(£)d£ as the probability that a randomly se- 
lected link has the length between £ and £+d£ and assume 
a power law distribution D(£) ~ I s . Since the network 
with a — has no length dependence and since the nodes 
are located in random positions in space with uniform 




0.0 0.2 0.4 / 0.6 0.8 




FIG. 2: (a) The probability density distribution D(£) of the 
link lengths I in d — 2 for five different a values: 1/2 (circle), 
(square), -1/2 (diamond), -1 (triangle up), -3/2 (triangle 
left) for networks of T = 2 12 (b) 5(a) varies linearly with a 
for a > —4 in d = 2 and for a > — 2 in d = 1. The saturation 
values of <5 are -3 and -2 for d — 2 (filled circle) and d = 1 
(opaque circle) respectively, (c) Scaling of the probability 
density for a = -oo and for T = 2 11 , 2 13 and 2 15 in d = 2. 

probabilities, D(£) for a = should depend only on the 
volume of the spherical shell between £ and £ + d£ and 
D(£) ~ £ d ~ 1 in the <i-dimensional space. For a^0, this 
distribution is modified by the factor £ a of Eqn. 3 giving 

D(£) - (4) 

which implies 5(a) = a + d—1. Therefore at a particular 
value of a = a c = 1 — d, 5(a) — and the distribution 
is uniform in any dimension. D(£) grows with £ for a > 
a c and decays with £ for a < a c . Growing, uniform 
and decaying distributions are shown in Fig. 2(a) for 
different values of a which confirm a c = —1 in d = 2. 
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FIG. 3: For networks with T — 2 12 , variations of (a) the 
probability P(l) of a node of degree 1 (b) average length £(a) 
of a link with a. Opaque circles for d — 1 and the filled circles 
are for d — 2. 

Similar calculations in d = 1 show uniform distribution 
for a c = 0. 

Since we distribute nodes within the unit square and 
link lengths are measured using the periodic boundary- 
conditions along both the x and y directions, the dis- 
tance between any two nodes can be at most £ = 2~ 1//2 . 
Consequently all orientations of links of lengths up to 1/2 
are equally likely. However links of lengths greater than 
1/2 have to be oriented more towards to the diagonal di- 
rections i.e., y — ±x lines to properly fit in and therefore 
their orientations are not equally likely. This anisotropic 
effect is observed in d=2 when D(£) decays at a faster 
rate for 1/2 < £ < £ D since isotropy of the orientation of 
these links is lost and therefore they are less probable. In 
contrast, we do not see such a region in one dimension in 
agreement with the theoretical analysis. 

However, 6(a) saturates at a minimum value S m below 
a certain value of a which is calculated exactly in the 
limit of a — > — oo. Since, in this limit a new node is 
always connected to its closest neighbour, the probability 
that the t + 1-th node has a link of length between £ and 
£ + d£ is 

D t (l)dl = a£ d - 1 t(l-b£ d ) t - 1 d£ (5) 

where b£ d is the volume of a hypersphere in d dimensions 
and a = bd. Therefore for a network evolved up to a time 
T, the link length distribution of the whole network is 

D(£)d£ = Y,J =1 D t {£)d£ = a£ d - 1 d£Zf =1 t(l - M*)*' 1 . (6) 

In the limit T — > oo this series converges for large £ to 
D(£) ~ £~ d ~ 1 giving S m = — (d+1). However, for small £, 
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FIG. 4: Average degree of the node (k(t,T)) at time t for 
a two dimensional network of duration T — 2 14 multiplied 
by t 1 /' 2 is plotted for different a. For a = 3,2,1,-1 plots 
are parallel to the t axis for large t. The curves deviate for 
a = -3/2, -2 and for -3. 

D(£) = c\£ d ~ x — c^t ~ x ignoring higher order corrections 
where ci ~ 0{t 2 ) and c-i ~ 0(t 3 ) for large t. This implies 
that D(£) must have a maximum at £ ~ T^ x l d for all 
d > 1 as verified below by the scaling analysis. For one 
dimension however no such maximum is expected and 
the power law decay starts right from the small values of 
£. Again since S varies linearly with a as 5(a) = a + d— 1 
in general, the minimum value 5 m is attained at a = —2d 
and remains same for smaller values of a. 

Our numerical findings nicely support these results: we 
check that 6 saturates nearly at S m = —3 for a < — 4 for 
d = 2 and 5 m = -2 for a < -2 in d = 1 (Fig. 2(b)). We 
look at the distribution in more detail in d = 2 for small 
values of t. On decreasing a from - la maximum appears 
at I w 0.01 (Fig. 2(c)) i.e., as £ increases from zero the 
distribution grows very rapidly as a power law, reaches a 
maximum and then decays. For a particular value of a, 
D(£) scales nicely with the duration T of growth as: 

D(£) ~ T x l' 1 G(lT x l' 1 ). (7) 

In Fig. 2(c) we plot this collapsed data which fits very 
well to the following form of the scaling function: 

Q(x) =a'x b '/(x 2 + c'f- (8) 

Values of all the constants a',6',c',rf' of this expression 
are dependent on a. No such maximum in D(£) at the 
small values of £ is observed in one dimension. 

We studied the cumulative degree distribution F(k) — 
/ fc °° P(k)dk ~ fc 1-7 and assume the following scaling be- 
haviour: 

F(k,T) ~r r 'T{k/Tt) (9) 

where the scaling function F(x) — ► x 1-7 when x « 1 
and T{x) — > constant for x >> 1. This implies that 
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7 = 1 + rj/(. For example for a = — 1 at d=2 a good 
data collapse is obtained for r) = 1 and £ = 1/2 giving 
7 = 3 and the same result is obtained for all values of 
a > — 1. However for a < — 1 we see a stretched exponen- 
tial variation P(k) ~ exp(—k^^) and '(/'(a) increases to 
1 i.e., to a pure exponential form at a — —oo. It therefore 
appears that the transition from the stretched exponen- 
tial to the scale-free behaviour is perhaps taking place 
at the specific value of a c = — 1 in d—2 and in general 
at a c = 1 — d in the d dimension. However a similar 
study in one dimension shows all indications that a c is 
very likely to be around -0.5, certainly greater than -1 
but seems to be smaller than zero. Therefore we conjec- 
ture that a c — 1 — d though our numerical analysis in 
one dimention does not fully confirm this prediction. At 
the opposite limit of a = oo i.e. when each node is con- 
nected to its farthest neighbour, the degree distribution 
is found to be exponentially decaying and it appears that 
it happens only at a — oo since even at a = 40 we found 
scaling of the distribution. 

We note a few more important properties of this net- 
work near a c . The first moment of the degree distribution 
(k) is exactly 2 since the sum of the degrees of all nodes is 
2T counting each link twice where as the number of nodes 
is T + 1. As the cut-off of the degree distribution varies 

as T c the (k) is defined as kP(k)dk/ P(k)dk. 
Assuming 7 > 2 the mean is (7 — l)/(7 — 2) in the large 
T limit. This gives 7 = 3 correctly. Since each node has 
exactly one outgoing link, the mean number of incoming 
links averaged over all nodes must be equal to 1. The 
fraction of nodes with nonzero incoming links should be 
l-P(l). In Fig. 3(a) we plot the variation of -P(l) with 
a. This varies from 0.4777 at a = —00 to very close 
to 2/3 for a > a c . Therefore the mean number of in- 
coming links per node (ki n (a)) averaged over all nodes 
with k in ^ is T/[(l - P(1))T]=3 for a > a c . This 
is checked numerically. We also keep track of the vari- 
ation of average link length (1(a)) as the characteristic 
distance of this system which is shown in Fig. 3(b). Since 
for a < a c smaller links are more probable, (i{a)) is very 
small and approaches zero as T- 1 ' 2 . For a > a c , {1(a)) 
grows with a. Finally the fraction of anisotropic links 
f(£ > 1/2) also grows from zero at a — a c . 

The average degree of a node which is introduced at 
the time t to a network which has grown up to a time 
T is denoted by (kit, T)) and follows a power law varia- 
tion with t at large t values as: (k(t,T)) ~ t 13 . For BA 
network, f3 — 1/2 is obtained [[fl| and it is connected 
to the degree distribution exponent 7 by the relation: 
[3("f — 1) = 1. We have calculated (k(t, T)) in two di- 
mensions for different values of a and in Fig. 4 we plot 
(kit, TWt 1 / 2 vs. t. on a linear scale. We ohservc that for 



(3 = 1/2 for this range of a values. 

To summarize, we studied a growing random network 
where the attachment probability 7r.; to a previous node 
i depends jointly on the degree of the node ki as well 
as on the cv-th power of the link- length I as iti ~ kil a . 
By tuning a we find that for a < a c , the degree distri- 
bution of the resulting network is stretched exponential 
whereas for a > a c the network is scale-free. We also 
observe that the link length distribution follows a power 
law: D(t) ~ £ s for the whole range of the parameter a in 
contrast to the Waxman's exponential distribution Jl2| . 
The exponent S grows linearly with a for a > —2d and 
saturates at — (d+1) for a < —2d. Our interesting obser- 
vation is when a = —d the network has the property of 
real internet network where the link length distribution 
varies inversely with the link length. 
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